Open quantum system approach to single-molecule spectroscopy 
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In this paper, single-molecule spectroscopy experiments based on continuous laser excitation are 
characterized through an open quantum system approach. The evolution of the fluorophore system 
follows from an effective Hamiltonian microscopic dynamic where its characteristic parameters, i.e., 
its electric dipole, transition frequency, and Rabi frequency, as well as the quantization of the back- 
ground electromagnetic field and their mutual interaction, are defined in an extended Hilbert space 
associated to the different configurational states of the local nano-environment. After tracing out 
the electromagnetic field and the configurational states, the fluorophore density matrix is written 
in terms of a Lindblad rate equation. Observables associated to the scattered laser field, like op- 
tical spectrum, intensity-intensity correlation, and photon-counting statistics, are obtained from a 
quantum-electrodynamic calculation also based on the effective microscopic dynamic. In contrast 
with stochastic models, this approach allows to describe in a unified way both the full quantum 
nature of the scattered laser field as well as the classical nature of the environment fluctuations. 
By analyzing different processes such as spectral diffusion, lifetime fluctuations, and light assisted 
processes, we exemplify the power of the present approach. 

PACS numbers: 42.50.Ct, 42.50.Ar, 03.65.Yz 



I. INTRODUCTION 

In spite that in recent years many different experi- 
mental techniques have been established, single molecule 
spectroscopy (SMS) [l|, 0, H|, that is, the study of single 
nano-objects is predominantly realized by purely optical 
methods, i.e., by measuring the far-electromagnetic field 
scattered by the system when it is subjected to laser ra- 
diation. 

In contrast to standard quantum optical systems [4|, 
H, H, 0, B Hi, where the reservoir is only defined by 
the background (free) quantized electromagnetic field, in 
SMS the supporting nano-environment is highly struc- 
tured. Its definition and specific properties depend 
on each experimental setup. In fact, the local envi- 
ronment felt by the flu orop hore may involve cryogenic 
single-molecules [f(], EH Il2l|. m olecules at room temper- 
ature, either immobilized [l3l [Llj or diffusing in solu- 
tion [I5| , the solid state matrix supporting a nanocrystal 

or single bio- 
Even more, the 
proper definition of the background electromagnetic field 
may be modified because the local surrounding of the sys- 
tem may develop fluctuations in its dielectric properties 

dill. 

Each specific environment leads to different underly- 
ing physical or chemical processes that in turn modify 
the emission properties of the fluorophore. In most of 
the situations, it can be model as a two-level optical 
transition. The central theoretical problem of SMS is 
to relate the scattered laser field statistics with the un- 
derlying nanoscopic environment dynamic. Due to its 
complexity, its dynamical influence is usually taken into 
account by adding random elements to the characteristic 
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parameters of the system 0, 0, For example, spec- 
tral diffusion processes are taken into account by adding 
classical noise fluctuations to the transition frequency of 
the system, while lifetime fluctuations may be associated 
to transitions between different conformational (chemi- 
cal or physical) states of the environment. Stochastic 
Bloch equations [H, [H, HH, [36[ , stochastic decay rates 
[37l HI] , modulated reaction models [2!| [3i| or stochas- 
tic reaction coordinates [4(| are some of the associated 
theoretical models. On the other hand, radiation pat- 
terns whose statistics depend on the external laser power 
[H HE H3, HO, E3, HI H SI are modeled by introduc- 
ing extra states coupled incoherently to the upper level 
of the system (46[. 

The previous approaches are well-accepted and stan- 
dard theoretical tools for modeling SMS experiments. 
They provide a solid basis for describing different ex- 
perimental observables, such as those associated to the 
photon counting statistics. Nevertheless, as the influ- 
ence of both the background electromagnetic field and 
the reservoir fluctuations are represented in a unified way 
by a classical noise, in general it is not clear how to cal- 
culate arbitrary observables associated to the scattered 
laser field. This limitation is evident when considering, 
for example, the optical or absorption spectrum, which 
is defined in terms of the scattered field correlations. 

The previous drawback or limitation could be sur- 
passed if one is able to describe SMS experiments on 
the basis of a full open quantum system approach. By 
an open quantum systems approach we mean: (i) the 
possibility of describing both the fluorophore system and 
the reservoir fluctuations through a density matrix op- 
erator, whose evolution, i.e., its master equation, can 
be written and adapted to each specific situation, (ii) 
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Describing the quantum nature of the scattered electro- 
magnetic field through operators and providing a closed 
solution for their evolution, in such a way that a simple 
procedure for calculating arbitrary field correlations is 
established. Then, the calculation of observables like the 
optical spectrum follows straightforwardly, (iii) Charac- 
terizing the photon counting process through a Mandel 
formula and establishing a manageable analytical tool for 
the explicit calculation of the photon counting probabil- 
ities. The main goal of this paper is to demonstrate that 
it is possible to build up such kind of powerful and gen- 
eral approach, which not only recovers the predictions of 
standard approaches, but also allows to characterize, in 
different situations, arbitrary observables associated to 
the scattered quantum electromagnetic field. 

The scarce application of an open quantum system ap- 
proach to SMS has a clear origin. Due to the complexity 
of the underlying nano-environment, a microscopic de- 
scription, from where to deduce the system density ma- 
trix evolution, is lacking. We overcome this difficulty by 
noting that the noise fluctuations induced by the reser- 
voir [J, 0, y] can be associated to a coarse grained rep- 
resentation of its complex structure [47| , allowing us to 
write microscopic interactions that take into account its 
leading dynamical effects and at the same time are an- 
alytically manageable. From the effective microscopic 
dynamic the system and reservoir fluctuations result de- 
scribed in terms of a Lindblad rate equation. The the- 
oretical validity of these equations for describing non- 
Markovian ope n q uantum system dynamics was estab- 
lished in Refs. [48[ El]- The possibility of establishing a 
quantum-electrodynamic treatment of SMS experiments 
also relies on the results of Ref. where the quan- 

tum regression hypothesis was analyzed for Lindblad rate 
master equations. 

We remark that recent author contributions antici- 
pated the possibility of establishing the powerful for- 
malism developed in this paper. In Refs. (5lj . it was 
demonstrated that an anomalous fluorescence blinking 
phenomenon can be dynamically induced by a complex 
environment whose action can be described by a di- 
rect sum of Markovian sub-reservoirs. In Ref. [52|, for 
the same situation, observables associated to the pho- 
ton counting process and a mapping with triplet blinking 
models [HHg] were characterized. In Ref. [53|, we shown 
that fluorescence blinking patterns whose statistics de- 
pend on the external laser power can be model through 
an underlying tripartite interaction. Both situations lead 
to quantum master equations that correspond to partic- 
ular cases of the present general formalism. Since in both 
cases the calculations relied on specific generalizations of 
a quantum jump approach [54[, there is not a recipe for 
calculating the scattered field correlations. The present 
quantum-electrodynamic treatment also fills up this gap. 

The paper is outlined as follows. In Sec. II, the effec- 
tive Hamiltonian microscopic dynamic that defines the 
full approach is introduced and the evolution of the flu- 
orophore density matrix is obtained. In Sec. Ill, observ- 



ables associated to the scattered laser field are derived 
from the full microscopic approach. In Sec. IV, the pho- 
ton counting statistics is characterized through a Mandel 
formula and a generating function approach (55l . |56| . In 
Sec. V, the formalism is exemplified by analyzing differ- 
ent kind of environment fluctuations. Processes like spec- 
tral diffusion, lifetime fluctuations, and light assisted pro- 
cesses are explicitly characterized through the scattered 
field observables. In Sec. VI we give the conclusions. 



II. EFFECTIVE MICROSCOPIC DESCRIPTION 
AND DENSITY MATRIX EVOLUTION 

The total Hilbert space TLr associated to a SMS ex- 
periment is defined by the external product Ht = Tis <8> 
Tijr ® He, where each contribution denotes respectively 
the Hilbert space of: the fluorophore system, the back- 
ground electromagnetic field, and the rest of the degrees 
of freedom that define the local nano-reservoir felt by 
the system. In general, it is impossible to know the 
total microscopic dynamic of the reservoir B. In or- 
der to bypass this task, following an argument presented 
by van Kampen [4?| . we notice that the complexity of 
the reservoir may admits a simpler general description. 
Since the nano-reservoir can only be indirectly observed 
through the fluorophore system, its Hilbert space struc- 
ture can not be resolved beyond the experimental reso- 
lution. Therefore, it is split as Hb = ©^"T^Br, where 
each subspace TLb r is defined by the set of all quantum 
states that lead to the same system dynamic [47[ ■ As the 
reservoir may be characterized by inordinately dense as 
well as by discrete manifolds of energy levels, some TLb r 
may be of finite dimension. 

Clearly, the subspaces {WgjjJ^f define the maxi- 
mal information about the reservoir Hilbert space struc- 
ture that can be achieved from a SMS experiment. 
Then, our approach consists in representing the reser- 
voir through a set of coarse grained "configurational 
macrostates"{|i?)}^™1 x , (R\R') = Srr', each one be- 
ing associated to each subspace, and writing the to- 
tal microscopic dynamic in the effective Hilbert space 
T-Cff = 7is <8 "Hp ® Ttu- The configurational Hilbert 
space Tiii is expanded by the (unknown) configurational 
macrostates. The approach is closed after defining both, 
the field quantization and its interaction with the system 
in the total effective Hilbert space, and the self-dynamics 
of the configurational states. 



A. Electromagnetic field quantization 

The classical Maxwell equations [57j 

1 $B 

V-D = 0, yxE = -- — , (la) 
V-B = 0, VxH = +-5, (lb) 
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provide the basis for the quantization of the electro- 
magnetic field in a material media [H, The macro- 
scopic field vectors are related by the constitutive rela- 
tions D = eE, and B = /iH, where e and /i are the 
macroscopic material constants and c is the light veloc- 
ity in free space. For an absorptionless dielectric, the 
Maxwell equations lead to the classical Hamiltonian field 



(2) 



where V is the volume of integration. 

The canonical field quantization follows from the ex- 
pansion of the electric and magnetic field in normal 
modes and by associating to each component a creation 
and annihilation photon operator. Here, the domain of 
quantization V is defined by the local surrounding of 
the system and not by the full volume of the support- 
ing media. In order to take into account the possible 
dependence of the local field on the configurational bath 
states, the macroscopic material constants associated to 
the volume V are written as operators in TLu , 



J2^r\r)(r\- (3) 



Then, depending on the environment configurational 
state {l-R)}, the local field is quantized in a media char- 
acterized by the (real) constants {sr,ii r }. Consistently, 
as the system have nanoscopic dimensions we assume 
that these constants, in the local surrounding of the 
fluorophore (V), do not have any spatial dependence. 
After standard calculations [5|, both the electric E = 
E(+) + E<-) and magnetic H = H<+) + H("> field op- 
erators can be written in terms of positive and negative 
frequency contributions, each one related to the other by 
an Hermitian conjugation operation. They read 

h (.R) I/O 

EW = iJ2\R) (^lE(ir) ^ k ' r a kA , (4a) 



k.A 
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H(+> =iJ2\R) (^IE(^) 4^ kr a kA . (4b) 



k.A 



Here, we assumed that inside the volume V the nor- 
mal modes can be well approximated by plane waves. 
As usual, the creation and annihilation operators sat- 
isfy the commutation relation [akA,a k / A >] = <5 kk '<5 AA ', 
where k and A denote the wave vector and polarization of 
the quantized mode respectively. In agreement with the 
transversability of the electromagnetic field, the polar- 
ization of the magnetic field reads e kA = (kxe kA ), where 
k = k/|k|. r denotes the position vector. The dispersion 
relation associated to each configurational state is 



,(R) _ 



(5) 



From now on, we take [i R — * 1. As is well known [5j 
this assumption is always valid in an optical regime. 



B. System and field Hamiltonians 

The fluorophore system is defined by a two-level op- 
tical transition with natural frequency to a- The upper 
and lower states are denoted as \b) and \a) respectively. 
We take into account that, depending on the configura- 
tional state \R) , the natural frequency uja may be shifted 
(a) 

a quantity 5ui A . These parameters define the spectral 
shift of the system associated to each state of the reser- 
voir Therefore, the system and field Hamiltonians 
are written as 

H s = fkJ (X z /2, Hjr = E kA ^kAa kA a k A, (6) 

where Hf follows from Eqs. (0} and @. a z is the z- 
Pauli matrix written in the base {|a) , \ b}}. The frequency 
operator of the system is defined by 



w =w A + ^Suj^'IR) (R\, 



(R) 



R 

while for the field it reads 

cd kA = 

R 



(7) 



(8) 



The frequencies are defined by the dispersion rela- 
tions Eq. ©. 



C. Dipole-field interaction 

The natural decay of the fluorophore is induced by the 
coupling of its electric dipole with the quantized electric 
field. Their interaction can be written as 



E r 



d A , 



(9) 



where is the (electric) dipole operator associated to 
the optical transition and E rA is the electric field op- 
erator at the position r A of the system. To take into 
account the different configurational states of the reser- 
voir, the standard definition of the dipole operator d A 
Q is generalized as 



(10) 



R,R' 



where {d^Ri} are vectors (assumed real) with units of 
electric dipole. and a are respectively the raising and 
lowering operators acting on the states {\a) , \b)}. The 
diagonal contributions d R ji define the dipole associated 
to each configurational state. The nondiagonal contribu- 
tions, dfmi with R R' , take into account the possibility 
of coupling different configurational states (those of finite 
dimension) through system transitions. 

The dipole-field interaction Eq. ^j, in a rotating wave 
approximation [1,(1,0,11, lrom E( l s - O and reads 



zdip = i_j ( K flfl' crt \R) (R'\ ° kA 

k,A 

R,R' 



K^,a\R') (R\, 



(11) 
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where the interaction parameters are defined by 



k.A 

^RR' 



1/2 



,ikrA 



(ekA • d RR > 



(12) 



D. Environment configurational fluctuations 

The previous definitions effectively take into account 
the influence of the different configurational states of the 
environment on the fluorophore and the background elec- 
tromagnetic field as well as into their mutual interaction. 
Now, in agreement with the van Kampen argument [47 1 , 
and consistently with real specific situations [J @, y|j 
where the properties of the nano-environment fluctuates 
in time, the configurational states (associated to dense 
bath manifolds) must be endowed with a mechanism able 
to induce transitions between them. We remark that the 
results developed in Refs. [5l|, HH, [53| do not take into 
account neither rely on this extra dynamical effect, which 
in the context of an open quantum system approach can 
only be recover with the present treatment. 

In order to maintain a full microscopic description of 
the effective dynamics, here the (incoherent) transitions 
\R) <-> \R'} are introduced through the Hamiltonian 



H 



w 



Hi 



uw- 



(13) 



Hu is the free Hamiltonian of hi, Hyy define an environ- 
ment (W) responsible for the transitions \R) <-> \R') , 
while Hjjyj defines their mutual interaction. We assume 
that the states {|-R)} are the eigenvalues of Hu with 
eigenvalues hu R , and W is defined by a continuous set of 
arbitrary bosonic normal modes 

H u = Y JR t^R.\R)(R\, //-v Y., ( 14 ) 

bj and bj are the creation and annihilation operators of 
VV respectively. The Hamiltonian Hyw reads 

H uw = hJ2 (Xrw \R) (R'\ bj+x&v \R') (R\ 6}). (15) 

j,R,R' 
R^R' 

Then, the transitions \R) <-> \R') are assisted by the cre- 
ation or destruction of bosonic excitations in each of the 
modes of frequency ujj. 



E. System density matrix evolution 

The previous analysis allow to define the unitary dy- 
namic, (d/dt)p T (t) — — (i/h)[Hf, Pr(t)], of the density 
matrix pq-{t) associated to the Hilbert space TC/ = 
Tis ® ® T~(-u ■ The total Hamiltonian read 



H T = H S + H' u 



H 



T 



Hdi 



pi 



(16) 



where each contribution follows from Eqs. (O, (fTl"]) and 
(Ti"3"]) . p-r(t) describe the statistical dynamical behavior 



of the system, the electromagnetic field, and the config- 
urational states. The joint dynamic of the fluorophore 
and the configurational states is encoded in the density 
matrix Pg U {t), which follows after tracing out the de- 
grees of freedom of the background electromagnetic field 
T and the reservoir W, i.e., p$u(t) = Tr^-vv[pr(*)]- The 
system density matrix pg(t) — T^u[Psu(^)]t can always 
be written as 



(17) 



where the base was used for taking the trace over 

hi. The auxiliary states p R (t) = (R\ p su (t) \R) , define 
the system dynamic "given" that the reservoir is in the 
configurational state \R) . The dynamic of the configura- 
tional states follows from the populations 



P R (t)=Tr s [p R (t)}, 



(18) 



which provide the probability that the reservoir is in the 
configurational state \R) at time t. Therefore, these ob- 
jects encode the statistical properties of the noise fluc- 
tuations introduced in standard stochastic approaches 
[E El- The expression Eq. (fT5)) follows straightfor- 
wardly from the definition P R (t) = (R\ p u (t) \R) , where 
Pu(t)=Trs[p S u(t)]. 

Both, the electromagnetic field T and the bath W are 
assumed to be Markovian reservoirs. Then, their influ- 
ence can be described through a standard Born-Markov 
approximation. After some algebra the evolution of 
each state Pn(i) can be written as a Lindblad rate equa- 
tion [mmi 

= ^[HR,P R (t)} J R ({D,p R (t)} + - J[p R (t)}) 



E 

R> 



R'RPr^+Yj 

Vrr'Pw 

(t) (19) 

R' 

lRR'J[pR>{t% 



where {■•■}+ denotes an anticonmutation operation. 
The Hamiltonian reads 



jj f^R 
Hr = — — cr z , 



uj r = uj a + 



(20) 



and the remaining system operators are defined by 

D = ctV/2, J[»]=<r»o*. (21) 

The rates associated to the reservoir W are given by 

4>rw = Kgw(ARR>)\x R R i '\ 2 - 

Here, gw(u) is the density of states of W and A RR > = 
wr — lu R i are the transition frequencies of hi. The rates 
associated to the field T read (7^ = 7^) 



Irr 



, = 27r W 



dkff^'>(k)|« 



k,A |2 
RR 



S(\k\c-u A ). (22) 
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Here, g^(k) is the density of states of the quantized 
electromagnetic field associated to each configurational 
state. By writing 



3^(k)dk = 



co 2 Ve 3 J 2 



(23) 



where dfl is the solid angle differential and the contri- 
bution e^ 2 follows from the dispersion relation Eq. (O, 
after a standard integration Q it follows 



1 My ,2 1/2 

4^3 I'M e R - 



4tt 3ftc 3 1 



1 Mu ,2 1/2 
\ClRR'\ e R , . 

(24) 

For simplicity, it was assumed that the dipole vector Arri 
[defined by Eq. (flU|) ] can be written as (1rr/ = dRRid, 
where cIrri is the modulus of the dipole vectors, and the 
unit vector d does not depend on the coefficient R and R' . 
While this simplification forbid us to analyze the angular 
dependence of the scattered radiation, the general case 
can be easily worked out from the present treatment. 

The previous expressions rely on the condition Arri <C 
u>a, i.e., the transition frequencies of hi are much smaller 
than the optical frequency uja- To simplify the analy- 
sis, in Eq. (JTUJ) we have discarded any shift Hamiltonian 
contribution induced by the microscopic dynamic. Fur- 
thermore, since the fiuorophore S is an optical transi- 
tion, we have assumed that the average number of ther- 
mal excitations of the electromagnetic field at the char- 
acteristic frequency uja are much smaller than one, i.e., 
n(ojA, T) <C 1, where n(uj, T) — [exp(fiw/fcsT) — 1] , ks 
being the Boltzmann constant and T the temperature as- 
sociated to the electromagnetic field T . This inequality 
allows to discard in Eq. (|19[) the contributions that leads 
to thermal excitations in S. 

When the fiuorophore is subjected to a resonant exter- 
nal laser field of frequency ujl, the system Hamiltonian 
becomes H R — > H R + H iascT , with 



(25) 



As before, the operators cfi and &' are the raising and 
lowering operators acting on the states {|a) , |6)}. The 
system-laser detuning is given by 



The Rabi frequency f2o reads 

n = Y,ttR\R) (R\ 



(26) 



(27) 



where each Qr measure the system-laser coupling for 
each configurational state. 

The Lindblad rate equation Eq. (JTUJ) is the central re- 
sult of this section. It effectively describes the action 
of the nanoscopic fluctuating environment over the fiuo- 
rophore. The scheme of Figure 1 symbolically represents 
all processes associated to this equation. The first line of 




FIG. 1: Scheme associated to the fiuorophore evolution 
Eq. ([19]) . see text. 



Eq. (fl9|) describe the self-system dynamic for each con- 
figurational state \R) . The decay rate {7^}, the natural 
transition frequency {ojr}, as well as the Rabi frequency 
{f2/j} may depend on the bath state. The second line 
describe transitions between the configurational states 
(with rates {4>r'r}), whose dynamic does not depend on 
the state of the fiuorophore (property symbolically rep- 
resented in Fig. 1 by the surrounding ellipse). Finally, 
the third line describes configurational transitions that 
are attempted by a transition between the upper and 
lower states of the system (rates {7^/^})- In Sec. IV, 
a detailed analysis of each contribution and its associ- 
ated "kinetic environment evolution" [i.e., the dynamics 
of P R {t), Eq. (|T8])] is presented. 

As neither the initial conditions or the dynam- 
ics [Eq. (fT9|) ] introduce any coherence between the 



bath macrostates, their density matrix can be writ- 
ten as f>u(t) — ^2 R PR(t)\R) (R\ ■ From the relation 
Tr u [OuPu(t)] = Tru[Ou(t)pu(p)], it follows that the dy- 
namics of operators On (t) acting on Tin is classical and 
dictated by the probabilities PR,(t). On the other hand, 
the system dynamic arise after "tracing out" [Eq. p?|) ] all 
internal transitions between the configurational states of 
the reservoir (see Fig. 1). Therefore, the evolution of the 
system density matrix must be highly non-Markovian. In 
fact, taking into account the results presented in Ref. [HI, 
it follows 



dp s (t) 
dt 



£HPs{t) 



drUt-T)p s (r). 



(28) 



Cr defines the system unitary dynamic. The equation 
that defines the superoperator L can be explicitly written 
in a Laplace domain in terms of the propagator associated 
to Eq. (HID (see Eqs. (58) and (62) in Ref. [H). 



III. SCATTERED FIELD OBSERVABLES 

In SMS experiments, the fiuorophore dynamic is indi- 
rectly probed by subjecting the system to laser radiation 
[Eq. (|2"5j) ] and measuring the scattered electromagnetic 
field. Therefore, while the master Eq. (TiT)|) completely 
characterizes the system dynamic, one is mainly inter- 
ested in observables associated to the scattered laser field. 
In general, these observables can be written as a function 
of the electric field. For example, the flux of energy Sjf 
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per unit ar ea p er unit time (module of the Pointyng vec- 
tor) reads [57| 

SAt) = ~|E(*) x H(t)t| = ^^|E(i)| 2 , (29) 

where the second equality follows from the relation be- 
tween E(t) and H(t) [see Eq. (gj)]. 

The time evolution of the electric field operator follows 
from a Heisenberg evolution with respect to the total 
Hamiltonian Eq. (16]), i.e., (d/dt)E(t) = (i/h)[H T , E(t)]. 
Taking into account Eq. (j4a)l . the time dependence of 
E(i) can be obtained from the evolution of the operators 
JR) - JR) 



'k,\ 



\R) (R\ . We get 



^■ a< kx(' t ) ■ (R) (R)u\ \^ k.A* / Qn \ 



rff 



where 



XR) 



R> 



and the coefficients 



arc defined by 



Eqs. (J5|) and (fT2|l respectively. Furthermore, the SU op- 
erator (Jrir is defined by 



gr,r = a \R') (R\ 



(31) 



Eq. p0|) relies on a set of approximations consistent with 
the effective representation of the reservoir. Since the 
dielectric constant of each configurational state is well 
defined, any non-diagonal non-linear coupling between 
the field modes is discarded. Terms arising from R' u 
[Eq. (fl~3)) ] are also disregarded because they only intro- 
duce a small modification to the natural (optical) fre- 
quency of each mode. 

The dynamics of (t) can be written as the addition 
of two contributions, each one associated to the homoge- 
neous and inhomogeneous terms in Eq. (|30p . Then, the 
electric field is written as 0, Q 



E(t) = E f (t)+E s (t). 



(32) 



The contribution Ef (t) defines the free evolution of the 
field. In terms of positive and negative frequency contri- 
butions, E f (t) = E { f + \t) + E[ -) (f), it is defined by 



e[ + V) = * y \r)(r\ 



fku 



RM.X 



2Ve R 



OkA(0). 

(33) 

The scattered field contribution E s (t), associated to the 
inhomogeneous term in Eq. (f30|) . after some algebra [j| 
can be written as E s (t) = F>i + \t) + E&'(t), where 



E s +) (*) = 4^:^ x f ) x f l H d RW a R , R {t - t/c r ), 



RR> 



with f = r/|r|, r = |r|, and cr — c/^/Fr. As before 
[Eq. (|24p], for simplicity we assumed that d R i R = cIr/rcI. 



The previous expression allows to obtain the electric 
field in terms of operators defined in the system and con- 
figurational Hilbert spaces. Consistently with the effec- 
tive representation of the reservoir, it is rewritten as 



RR' 



(35) 



where each contribution reads 



E 



(+)(*) 

s(R\R'] 



4^[(d x f) x f] d R R,a R ,{t - r/c R ). (36) 

Since the background electromagnetic field [Eq. (|4a[)] 
does not involve any coherence between the configu- 
rational states, the operator [\R') (R\](t) appearing in 
Eq. ([33)1 must be read as follows. It labels all contri- 
butions to Eg (t) that at time t are attempted by the 
configurational transition \R) — > \R) , (R ^ R'). Then, 

the "conditional operator" E^S,, defines the electric 
field restricted to the condition that at time t the reser- 
voir is in the configurational state \R') and change to the 
state \R) . Similarly, the diagonal contributions E^fm 
define the electric field "given'" that at time t the reser- 
voir is in the configurational state \R) . Consistently, the 
conditional operator cr R (t) gives the evolution of the sys- 
tem operator a(t) under the same condition. 

Any scattered field observable must be written in terms 



of the conditional system operators E^^L 
characterize the field correlations. 



Below, we 



A. Correlations 

When the scattered field is measured with photo- 
electric detectors, the usual observables can be writ- 
ten in terms of two time (normal order) correlations 

as, sal 



d(r) = lim : E$ -) (t)E$ +) (t + t) :, (37) 

t— *0O 



as well as 
C 2 (r 



lim : E< _) (t)Er> (t + r)Er ; (t + rjE^' (t) 

t — >oo 



*(+)/ 



(38) 

The overbar denotes quantum average with respect to the 
total initial density matrix. The symbol : (• • • ) : takes 
into account the right interpretation of Eq. ([34]) and de- 
notes a summation over all internal configurational paths 
defined through the conditional contributions E^i^n, 
Eq. (|36| . Then, the correlations are explicitly written as 

\ ' £(-)(*) 



C 1 (r)= lim J2 



L^i ^ S (-R 2 |fli) s(-R.4|-R 3 )' 



(39) 



R,!,R 2 R 3 ,Ri 



and 



C 2 (t) = lim Y E 



£.(-)(*) p(-)(i+-r)^(+)(t+r) £,(+)(*) 

' (R 2 \Ri)' 



t^oo^ '- l s(R 2 \Ri) rj s(R i \R 3 ) rj s (R i \R 3 ) 



Ri,R-2 
R 3 ,Ra 



(40) 
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where E 



(±)(t) 



E 



(±)(«) ..1/4 



J s(^|fl') - " s (ii|fl') £ K' • This definition guaran- 
tees that both correlations can be related to observ- 
ables defined in terms of energy (photon) fluxes [see 
Eq. (|29|) ] . It also takes into account that before the tran- 
sition \R') — > \R) , the dielectric constant of the bath is 
£r>- 

By using Eq. (f36|) and the rate expressions Eq. (|24)l . 
the first order correlation can be written as 



Ci(r) 



/(r 



Ancr 2 ' 



where for brevity we define /(f 

1r = Ir + ^Ir'R- 



(41) 

[(d x f ) x f] 2 , and 



(42) 



Eq. (|4"Tj) can be considered as a natural generalization 
of the expression corresponding to the Markovian case 
[58| . Similarly, the correlation of the (conditional) rais- 
ing and lowering operators can also be obtained from the 
system density matrix evolution [Eq. (|19p ] after invok- 
ing to a quantum regression theorem. For Lindblad rate 
equations it reads [50J 

Ox (t)A(t + r)0 2 {t) = J2 Trs{A(e Tl ) RR , [0 2 p R , (t)Oi]}, 

RR' 

(43) 

where 0\, O2 and A are arbitrary system operators. Each 
contribution [indexed by R, R'] , defines a conditional av- 
erage of the involved operators: at time t, the configura- 
tional bath state is \R'} , while at time (t + r), it is \R) . 
C denotes the generator of the Lindblad rate equation, 



e tC ) RR >p R ,(0). 



(44) 



From Eqs. (gSj) and Eq. {41]) it follows 



Ci(r)=/(f) 



Aircr 2 

RR' 



s/MR~^sW{e rC )Rw\p^}}, 

(45) 

= lim^oo p R [t). By using the same calcula- 
tions steps, the second order correlation reads 



C 2 (r) = /(f) 



Alter 2 J ^ 

RR' 



7 fl Tr s {<7V(e T£ W (46) 



iR'VPR'^ 



J21R'R"< j PR"< j 
R" 



The expressions Eqs. (|4"5"]) and (|4"6"j) are the central results 
of this section. They allow to characterize observables 
such as the spectrum of the radiated field as well as the 
intensity-intensity correlation. 



B. Spectrum 

The spectral intensity radiation (in units of energy 
ftwA) per unit of solid angle [J, Q is defined by the di- 



mensionless expression 



5 M = or 

Z7T 



+ 00 



dre 



where ldl is the frequency of the laser excitation, 
[Eq. (|25p]. As usually lim^oo Ci(r) ^ 0, the spectrum 
can be split in a coherent and incoherent contributions 

S(uj) = f(v)[S coh (Lo) + (2Tr)- 1 S mc (uj)}. (48) 

S C oh{w) consists in a Dirac delta term 



S coh (u) = sfJ h 6(ujL) 



(49) 



that measures the scattered radiation emitted at the fre- 
quency of the laser excitation. From Eq. (145 p . it follows 

S ( l = £ (b\ PR \a) (a\ p% \b) . (50) 

RR' 

The incoherent contribution can be written as 

Sinc(w) = Ci(u)\ + C 1 {u)\ (51) 

where C\ (u) = dre~ UT Ci (t) is the Laplace transform 
of the function C*i(r) = [ftw A /(f)/47rcr 2 ]- 1 Ci(r). 

In a Markovian limit, i.e., when the configurational 
space is unidimensional, from Eq. (|5ip it is possible to 
write the spectrum as an addition of three Lorentzian 
functions (Mollow triplet) whose widths and heights de- 
pend on the natural decay of the system and the Rabi 
frequency @, [E H • I n tne general non-Markovian case, 
the spectrum also has a strong dependence on the param- 
eters that define the environment fluctuations (see next 
sections). In spite of these dissimilarities, in both cases 
the spectrum is mainly related to the dynamic behavior 
of the system coherences [see Eq. ([50]) ]. 



C. Intensity-Intensity correlation 

The normalized intensity-intensity correlation Q reads 



:I(t + r)I(t) : C 2 (r) 
g 2 {T) = hm ==s = — (52) 

t^oo ^ \Cl(0)\ 2 

where I(t) = ~Ei~\t)~Ei + \t) is the intensity operator. 
From Eqs. {15]) and {15]) it follows 

= 7A2 E^H ^ ^ Tt ) RR .[at,] \b) , (53) 

(1st) 



RR' 



where 



*S = [lR> (b\PR'\b)+Y,7*n" WPk"^ }l«X«l > (54) 
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and the normalization constant reads 



E^ WPnW 



(55) 



R 



As in the Markovian case [j| , Eq. (|53|) corresponds to the 
probability density of detecting one photon in the sta- 
tionary regime (lim^oo) and a second one in the interval 
(t, r + dr). The factor a^, takes into account all possible 
emission paths that leave the system in the ground state 
and the reservoir in the configurational state \R') . The 
sum over the index R takes into account all photon emis- 
sions that happen in the interval (r,r + dr) and leave 
the reservoir in the state \R) . The normalization factor 
I s t defines the average stationary intensity emitted by the 
fluorophore. 



IV. PHOTON COUNTING STATISTICS: A 
GENERATING FUNCTION APPROACH 

In most of the SMS experiments the scattered laser 
radiation is measured with photon detectors. Then, the 
photon counting statistics is also an usual observable. As 
in standard fluorescent systems, the probability P n (t) of 
detectingn photons up to time t follows from a Mandel 
formula 0, Q- Here, it is generalized as 



Pn{t) = ~ 





n 




[ dt'l(t') 

'-Jo 


exp 


- f dt'I(t') 
Jo - 1 



(56) 



where the normalized intensity operator I(t) = 
{hLU A /4:ircr%)- 1 ~E,i~ ) {t)~E,i +) (t) has units of photon flux. 
rd denotes the distance between the fluorophore and the 
detector. : (• • • ) : denotes both an usual (normal) time 
ordering and a summation over all internal configura- 
tional paths, whose definition follows straightforwardly 
from Eqs. and ([iO ]) . 

While Eq. (fB"6"]) allows to characterize the probabilities 
P n (t), a simpler technique to calculate these objects is 
provided by a generating function approach [55|, 56] . This 
very well known technique was also used in the context of 
SMS when dealing with stochastic Bloch equations [H, 
[HI, [H, [36| and related approximations. In contrast, here 
we formulate the generating function approach [55|,[5(| on 
the basis of one of the central results of this contribution, 
i.e., Eq. (fTO|) . Added to its broad generality, the present 
formulation avoid the use of any stochastic calculus. 

By writing the system density matrix as 



— ^n— 



(57) 



where each state (t) corresponds to the system state 
conditioned to n— photon detection events p, |54|, the 
probability of counting n-photons up to time t reads 



P„(t) = TV s [pW(t)]. 



(58) 



A "generating operator" [55| is defined by 

* — 'n— 



(59) 



where s is an extra real parameter. This operator also 
encodes the system dynamic, Pg(t) = G(t, s)\ a= \. The 
conditional states p^ n \t) can be decomposed into the 
contributions associated to each configurational state of 
the reservoir, leading to the expression 



G{t,s) 



£ 



n=0 



( n )/v> - 



G R (t,a). (60) 



Each matrix p^ (t) defines the state of the system under 
the condition that, at time t, n— photon detection events 
happened and the configurational state of the environ- 
ment is \R) . Consistently, each contribution Gn(t, s) de- 
fines the (conditional) generating operator "given" that 
the reservoir is in the configurational state \R) . Its evo- 
lution, from Eq. (|19p . can be written as 

dG R (t,s) 



dt 



— \H R , G R (t, a)] - i R {D, G R (t, s)}+ 



Y,4> R ' R G R {t,s) (61) 

W R' 

+s lR J[G R (t, a)] + s E 1rr?J[Gr> (t, «)], 

R' 

where 7 fl is defined by Eq. (f42|) . Notice that the param- 
eter s is introduced in all terms (third line) associated 
to a photon detection event, i.e., those proportional to 
oW |,[5|. 

In the context of SMS [H S, S ^ , the matrix el- 
ements of G(t,s), in an interaction representation with 
respect to HlolO'z, are usually denoted as 



U{t,s) = 


\{G ab e 




V(t,s) = 


±(G ab e 


—iu>Lt 


W(t,a) = 




Gaa) ) 


y(t,s) = 


2^ Ghh + 





b a e 



G 



b a e 



(62a) 
(62b) 
(62c) 
(62d) 



[G pq = (p\ G(t, s) \q}], while their evolution is called "gen- 
eralized optical Bloch equation." From Eq. (|6"U|) . it is pos- 
sible to write each matrix element as a sum over the pa- 
rameter R. In Appendix A, we provide the evolution of 
each component [Eq. (|A1[> ]. 

After getting the matrix elements of the generating 
operator, Eq. (|62|) . the photon counting process can be 
characterized in a standard way [47j . From the definition 
of the generating operator, Eq. (f59"|) , the photon counting 
probabilities Eq. ([55]) follows as 



(63) 



s=0 
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The function y(t,s) also allows to calculate all factorial 
moments 

oo 

N (k) {t) = ^ n(n - 1) • • • (n - k + l)P n (t), (64) 

n=0 

in terms of its derivatives 



^ fc) (*) = 2 5? y(*,*) 



(65) 



Furthermore, the first two moments of the photon count- 
ing process, N k (t) = J2n=a nkp n(t), {k = 1,2), read 



AT(t) = 2^-y{t,s) 
os 

^(*) = 2 £^' s ) 



2^-y(t,«) 
as 



(66a) 
(66b) 



Both moments encode important information about the 
scattered radiation. The line shape of the fluorophorc 
system is defined by [l| 



IK) = lim -N(t), 
t— *oo at 

while the Mandel factor is defined as 

A^( t) -lV 3 (t) _ 
Wj " AT(t) 



(67) 



(68) 



As is well known 0, M, @, 0] , it allows to determining the 
sub- or super-Poissonian character of the photon count- 
ing process. In Appendix B we show the consistency 
between the generating function approach, the Mandel 
formula Eq. (|56p , and the results obtained in the previous 
section. In Appendix C, it is shown that the stationary 
Mandel factor 



Q st = lim Q(t), 

t — >oc 



(69) 



can be obtained in an exact analytical way after solving 
the evolution Eq. (fBTj) in the Laplace domain. 



V. EXAMPLES 

In this section, different processes covered by Eq. (fl9|) 
are analyzed. The examples are classified in accordance 
with the evolution of the environment fluctuations, i.e., 
the evolution of the configurational populations Pn(t), 
Eq. (fTSf . In each case, observables such as the spec- 
trum Eq. (|4"T|. intensity-intensity correlation, Eq. (|55|) . 
line shape, Eq. (f6"7) , and Mandel factor, Eq. can 
be calculated in an exact analytical way. In fact, in a 
Laplace domain, both the evolution of the auxiliary den- 
sity states, Eq. (TT9"]) , and the evolution of the auxiliary 
generating operators, Eq. (|6ip . become algebraic linear 
equations. 



A. Self environment fluctuations 

When the nano-reservoir only has dense manifold of 
states, the fluctuations between the configurational states 
must be governed by a classical master equation [47l |. 
This case is covered by taking 



dt 



-Wr, p R (t)] - 1r ({D, p R (t)}+ - J[p R (t)]) 



^r'rPrW) + Y <Prr>Pr> (*)) 

R' R' 



(70) 



which arises from Eq. (|19| under the condition d^< = 0, 
R 7^ R', implying the vanishing of the rates {7^'} 
[Eq. (f2g))]. From Eq. K7DJ), the evolution of the config- 
urational populations, Eq. (fl8|) . reads 



<pRR'PR'(t)- ( 71 ) 



In the context of SMS, this evolution defines the "ki- 
netic dynamic" of the environment. Consistently, this 
equation, and in consequence the transitions between the 
configurational states {|-R)}^™i x , does not depend on the 
state of the system. Notice that in the effective micro- 
scopic representation, the configurational transitions are 
induced by the reservoir W, Eq. (fT3)) . 

While the fluctuations between the configurational 
states are defined by the classical evolution Eq. ([7T|) , the 
quantum dynamic of the system, for each bath state, is 
defined by the diagonal contribution in Eq. (|70|) . De- 
pending on its properties, different processes are recov- 
ered. 



1. Spectral fluctuations 

A nano-environment consisting on molecules at very 
low temperatures may induce random fluctuations in the 
natural frequency of a fluorophore [l|. This situation 
can be described by the present formalism by assum- 
ing that the natural frequency of the system, Eq. ([^U]) . 
is the unique parameter that depends on the configura- 
tional bath states. Consequently, each configurational 
bath state has associated a different system transition 
frequency, u> R = uja + Slur. Equivalently, each configu- 
rational state "induces" a spectral shift defined by the 
parameter 8ojr. The population Pr{€) gives the proba- 
bility that, at time t, the natural frequency of the system 
is lor. The coefficients associated to the Rabi frequency, 
Eq. (|27[) . must be chosen all the same. Similarly, the de- 
cay rates associated to each bath state, Eq. (l24l) , are also 
all the same. Therefore, the system is characterized by 
one single Rabi frequency and a natural decay rate, i.e., 



R 



1r = 7- 



(72) 



In the following analysis, the configurational space is 
assumed two dimensional, R = 1,2. Furthermore, the 
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FIG. 2: Optical spectra [Eq. (|5ip ] for a fluorophore system 
subject to spectral fluctuations in the regime (/> <ti Suj 
7 ~ fi. (a) The parameters, in units of the decay rate 7, are 
Q. = j/y2, 5ui — 7/IO, and <f> = 7/125. (b) The parameters 
are Q — 57, 8ui = 7/IO, and cj> — 7/5OO. The inset shows the 
narrow spectral peak. In both cases, the laser frequency is in 
resonance with the system, lul — u>a- 

problem is restricted to a symmetrical one. The charac- 
teristic parameters read 

8u>^ — +5u>, 4> 2 i — 4>i (73a) 

Sujf = -Suj, 12 = 0. (73b) 

With these assumptions, from the results of Appendix C, 
it is possible to characterize both the line shape, Eq. (ft!?)) , 
and the stationary Mandel factor, Eq. (|69|) . in an exact 
analytical way. The final expressions recover the results 
obtained by He and Barkai in Ref. [36(, where the nat- 
ural frequency of the fluorophore is modeled through a 
stochastic two-state process. 

Now, going beyond the possibilities of any previous 
approach, in the following figures we characterize the 
spectrum of the scattered radiation [Eq. (|5T|) ]. Its ex- 
act analytical expression is not provided due to its ex- 
tension. In all cases, the system-laser detuning [Eq. (f2T))) ] 
is zero, i.e., ujl — ^>a- When the characteristic time be- 
tween configurational transitions, (ff~ , is the small time 
scale of the problem, the spectrum shape strongly de- 
pends on the relation between the natural decay 7, the 



Rabi frequency fl, and the spectral shift Suj. In Fig. 2, 
the condition (f> <C Suj <C 7 ~ f2 is satisfied. Then, the 
spectral fluctuations induced by the bath slightly mod- 
ify quantum dynamic of the system. Nevertheless, the 
spectrum shape presents strong deviations with respect 
to the Markovian case. 

In Fig. 2a, as £1 < 7, the spectrum only consists in 
one Lorentzian component (Rayleigh peak) [54[. In con- 
trast to standard (Markovian) fluorescent systems, here 
a central narrow peak is clearly visible. It must not be 
confused with the coherent Dirac delta contribution that 
appears in both Markovian and non-Markovian systems 
[Eq. 1(15])]. In fact, as in three level systems its origin 
can be related to the dynamics of the system coherences 
[Hi^ ]. The spectrum develops a narrow peak whenever 
the dynamic of the coherences slowly fluctuates between 
two different regimes. The width of the narrow peak is 
given by the addition of the constant rates that define 
the fluctuations between the different dynamical regimes 
[54| . Here, as <fi <C Suj, each dynamical regime is defined 
by the spectral shifts ±5uj associated to each configu- 
rational bath state. The "blinking" between these two 
regimes is governed by the rate 4>. Therefore, the width 
of the narrow peak is 2<fi, which allows us to read a cen- 
tral property of the environment fluctuations from the 
optical spectra. 

In Markovian fluorescent systems, when fl > 7 the 
laser strongly couples the population and coherences of 
the system, leading to oscillations that induce two extra 
peaks in the optical spectra (Mollow triplet) Q. Here, 
this phenomenon is shown in Fig. 2b, where the param- 
eters satisfy <fi <C Suj <C 7 < ft. As the Rabi frequency 
n is much larger than the spectral shifts Suj, the peaks 
appear at ±J7. Consistently with the parameter values, 
the extra narrow peak (inset) is also associated with the 
coherences blinking-like evolution. 

In Fig. 3a, the condition <p <C 7 — Vt <C Suj is satisfied. 
Therefore, the quantum system dynamic is mainly gov- 
erned by the spectral shifts induced by the environment. 
Consistently, the spectrum develops two peaks at ± Suj. 
In contrast with the Mollow triplet, which is associated 
to the Rabi oscillations, here the central Rayleigh peak 
(lul = loa) is almost absent. It only appears a narrow 
central peak (inset) related to the coherence dynamics. 

By maintaining the value of all parameters, in Fig. 3b, 
the laser intensity was increased such that « 7 C (] ~ 
Suj. In this situation, there is a competition between the 
Rabi oscillations and the spectral shifts induced by the 
laser and the reservoir respectively. Notice that the ex- 
tra peaks do not appear at ±f2 [Fig. 2b] neither at ±Su> 
[Fig. 3a]. Furthermore, in contrast to standard fluores- 
cent systems, the height of the central Rayleigh peak is 
smaller than the height of the lateral peaks. The inset 
shows the narrow central peak. Since the parameters 
of the environmental fluctuations are the same than in 
Fig. 3a, the central narrow peaks have the same widths, 
i.e., 2(p. 

A broad kind of spectrum behavior may arise when 
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FIG. 3: Optical spectra [Eq. (I5ip ] for a fluorophore system 
subject to spectral fluctuations in the regimes ^ < 7 ~ S] < 
5u> (a) and (f> <C 7 <C Q ~ <5oj (b). In (a) the parameters, 
in units of the decay rate 7, are f2 = j/\ // 2, <5w = 57, and 
4> — 7/(4 x 10 4 ). In (b) the parameters are Q. = 57, <5u = 57, 
and <j> = 7/(4 x 10 4 ). In both cases, the laser frequency is in 
resonance with the system, ujl = Ui. The insets corresponds 
to the central narrow peaks. 




FIG. 4: Optical spectra [Eq. (|51[) ] f° r a fluorophore system 
subject to spectral fluctuations. In (a) the parameters, in 
units of the decay rate 7, are Q = j/y2, Su> — 57, and <j> = 
IO7 (dashed line), cf> — 5O7 (dotted line), (f> = 1257 (solid line). 
In (b) the parameters are Q = 7, 5u) = 37, and <j> = 0.257 
(dashed line), = 7 (dotted line), (f> — 47 (solid line). In both 
cases, the laser frequency is in resonance with the system, 

LU L = UJA- 



the environment fluctuation time 4>~ is not the small 
time scale of the problem. In this regime, the coherences 
dynamic loose its dichotomic character and consequently 
the spectrum docs not develop any central narrow peak 
[59l ] . In Fig. 4a the reservoir fluctuations are faster than 
the system dynamic, 7 ~ ~ 5u> <C 4>. The spectrum de- 
velops the well-known effect of motional narrowing, i.e., 
by increasing the rate (f> the optical spectrum becomes 
narrow. In Fig. 4b the rate <f> is increased around the 
intermediate regime 7 ~ (f>. As the parameters of the sys- 
tem and the environment fluctuations are of the same 
order, a small variation of (f> lead to strong changes in 
the spectral shape. This property can also be found in 
the photon counting statistics (see Fig. 5 in Ref. [36j]). 



2. Lifetime fluctuations 

Chemical or physical conformational changes of the 
supporting nano-environment may modify the lifetime of 
the fluorophore [H, 0, This situation is covered by 



taking (j RR , = 0) 

n R - f>, Slur = 0. (74) 

Then, each configurational state \R) has associated the 
natural decay j R . The classical evolution Eq. (f7Tj) de- 
scribe the transitions, with rates 4>r>r, between the dif- 
ferent decay rates felt by the fluorophore. 



3. Molecules diffusing in solution 

The present formalism may also be applied for describ- 
ing the statistical properties of radiation patterns pro- 
duced by molecules diffusing in a solution [2|, HH • If the 
intensity of the laser field varies appreciably along the 
diffusing space, the Rabi frequency (which depends lin- 
early on the external laser power) must be considered 
dependent on position. Then, the configurational space 
can be associated to the physical space where the diffu- 
sion process happens. This situation is covered by taking 



12 



(7 



RR' 



0) 



1r = 7, 



Suji 



0. 



(75) 



Each coefficient il/j [Eq. ([27]) ] measures the laser intensity 
at each position labeled by the index R. The evolution 
Eq. (|7T|) describes the molecule diffusion. For an homo- 
geneous coupling between first neighbors, <fi RR ± 1 = <f>, 
a standard diffusion process with diffusion coefficient 
Dq = Ax 2 4>, is recovered. Aa; is the discretization step 
in physical space. 



B. Light assisted environment fluctuations 

When the local nano-reservoir surrounding the system 
only consists in a few degrees of freedom, its Hilbert space 
is defined by a discrete (finite) set of states. Therefore, 
each one may represents a different configurational state. 
While there exist different Lindblad rate equations that 
may model this situation (see end of this section), here 
we consider the case 



d PR(t) 
dt 



— [H R ,p R {t)]~ lR {{D,p R {t)} + -J[p R (t)]) 



This evolution follows from Eq. ([HI)) after taking 
{4> RR /} = 0, condition consistent with the finite structure 
of the reservoir, i.e., its dynamic is unable to induces self 
(internal) fluctuations [43]. In the effective microscopic 
representation, it follows from the interaction Eq. (|16[) by 
uncoupling the configurational space from the reservoir 
W, i.e., H uw = [Eq. (HID]. 

Taking into account the interaction Hamiltonian 
Eq. (fTTj) . we deduce that transitions between the differ- 
ent configurational bath states may only happen when 
they are attempted by a photon emission process, i.e., 
a transition between the upper and lower system states 
(see Fig. 1). Therefore, in contrast with the previous 
case, here the dynamic of the configurational populations 
[Eq. (fl~8|) ] strongly depends on the state of the system. 
From Eq. ([75]) we get 



+I> BIi .<(<), (77) 

R' R' 

where P R b \t) = (b\ p R (t) \b) . The evolution of P R \t) also 
follows from Eq. (|76p , which in turn involves the remain- 
ing matrix elements of all auxiliary states. Therefore, in 
general it is not possible to write a simple equation for the 
evolution of the configurational populations. Neverthe- 
less, in the limits {j RR ,} < {7^} and {j RR ,} > {7^}, 
an evolution similar to Eq. ([7T|) can be obtained. We 
assume that 



n, 



n. 



Suj r = 0, 



(78) 



implying that the fluorophore, independently of the con- 
figurational state, always feels the same laser intensity 
il, and maintains its transition frequency uja- These con- 
ditions allow to model a class of light assisted processes 
consistent with different experimental situations [53| . 

Under the condition {7^/} <C {"f R }, from Eq. ([76]) 
it is simple to realize that before happening a configura- 
tional transition, the fluorophore may emit a large num- 
ber of photons. Thus, "while" the bath remains in the 
configurational state \R) , the fluorophore can be well ap- 
proximated by a Markovian fluorescent system with de- 
cay rate 7^ [Eq. (|4"2"| ] , the radiation pattern being char- 
acterized by the average intensity 



Ir = 



IrP? 



1R 



2fi 2 + W 



(79) 



where S = uj^—oja is the system-laser detuning, Eq. ([26 
Consequently, it is possible to approximate 



(80) 



The first term represents the photon flux produced by 
the system while the bath remains in the configurational 
state \R) . This quantity is approximated by the intensity 
I R multiplied by the probability of being in the configura- 
tional state \R) , i.e., Pn(t). Similarly, a straightforward 
interpretation of this approximation follows by noticing 
that IrI^ r gives the stationary upper population of a 
Markovian optical transition with characteristic decay 
rate 7^. Then, Eq. ([50)) implies that, "given" that the 
reservoir remains in the state \R) , the upper population 

P R (t) quickly reaches the stationary value I R /~i R . 

By introducing Eq. (JSOj) in the evolution Eq. (J77J, it 
follows the classical evolution 



^RR'PR'{t), (81) 



where the transition rates read 



Tr'r = 7/?/ 



n 2 



2n 2 



4S 2 



(82) 



These expressions generalize the results presented in 
Ref. [HI]. The associated radiation pattern develops a 
blinking phenomenon, i.e., the intensity of the scattered 
radiation randomly changes between the set of values 
{Ir} associated to each configurational state. These 
changes are governed by the classical evolution Eq. (|81|) . 
The dependence of the transition rates T R i R on both the 
laser intensity and the system-laser detuning [compare 
with Eq. (|71[)]. gives the light assisted character of the 
blinking phenomenon. In contrast with triplet blinking 
models [l|, HH , here all rates that govern the transitions 
between the configurational states (intensity states) de- 
pend on both, the laser intensity and the system laser 
detuning, the first property being consistent with differ- 
ent experimental results [53| . 
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FIG. 5: Optical spectra [Eq. (|51|l ] (a) and intensity-intensity 
correlation [Eq. (|53p ] (b) for a light assisted blinking radiation 
pattern. The configurational space is two dimensional, R — 
1, 2. The parameter values, in units of the Rabi frequency SI, 
are j 1 = Q, 7 2 = WQ, j 21 = 0. 00150, j 12 = 0.02O. The laser 
is in resonance with the system, ujl — ua- In (a) the inset 
shows the central narrow peak. 



When {jrri} S> {Jr}, in the intermediate time be- 
tween two consecutive configurational transitions, the 
system is only able to emits a few number of photons. 
Therefore, the blinking phenomenon is lost because it is 
not possible to associate an intensity regime [with value 
Eq. (|79p] to each configurational state. In spite of this 
fact, under the replacement 7^ — > (J2r>> 1r"r) the ap- 
proximation Eq. (|80|) is still valid. Since {7^'} {7_r}, 
it follows j R ~ (J2r" 1r"r)i implying that both Eq. (j£Tj) 
and l|82p remain valid when {7^/} 3> {7_r}- 

In the following figures, the light assisted blinking phe- 
nomena described in Ref. [53[ is characterized through 
the scattered field observables. A two dimensional con- 
figurational space is assumed (R = 1,2) and also the 
condition {^rr>} <C {7^} is satisfied. Fig. 5a shows the 
optical spectrum [Eq. (|5Tj) ]. The Rayleigh central peak 
is endowed with a narrow peak (inset). As in the previ- 
ous analysis, it origin relics on the blinking like behavior 
of the system coherences, effect associated to Eq. (j5Tj). 
The width of the narrow peak (59)], T12 + I^i, depends 
on both the system-laser detuning and the laser intensity 
[see Eq. (1521) ]. 

Fig. 5b shows the normalized intensity-intensity cor- 




(co L -co A )/Q 



FIG. 6: Stationary Mandel factor [Eq. (|69[l ] as a function of 
the system-laser detuning. For the light assisted process (full 
line) the parameters are the same than in Fig. 5. The in- 
set shows the associated line shape [Eq. (|67[l ]. For the self 
fluctuating environment (dotted line) the parameters are de- 
fined from the mapping Eq. (|83p . The dashed line corre- 
sponds to the scaling Eq. (|85|l with So — Q, fl — 0.250, and 
7 12 = 0.0070. 



relation [Eq. ([53]) ]. For short times, this function satis- 
fies < #2 (t) < 1, property associated to photon an- 
tibunching. At intermediates times the correlation sat- 
isfies <?2 (t) > 1, indicating photon bunching, which in 
turn is a direct manifestation of the environment fluctu- 
ations. Furthermore, 52(7") is almost constant during a 
long period of time. This property is related with the 
telegraphic nature (high and low intensity states) of the 
scattered radiation field 0, [54[ • 

The central narrow peak in the spectra as well as 
the photon bunching- antibunching transition reflected by 
the intensity-intensity correlation (Fig. 5) are associated 
to the blinking property of the scattered field intensity. 
Then, it is expected that, independently of the under- 
lying dynamic, these features also arise whenever the 
radiated field fluctuates between two different intensity 
regimes. For example, the blinking phenomenon defined 
by Eq. (fgTj) and (f52"| . also can arises from a fluorophorc 
whose environment, independently of the system state 
[Eq. I|7ip]. randomly changes its natural decay between 
two different values, i.e., Eq. (|70[) with the parameters 
values Eq. I|74[) . Consistently, we have checked that un- 
der the mapping 



7i 

72 



7i 

72 



(7i 
(7 2 



■72iJ 
712) 



6 21 =r 21 , (83a) 
6 12 =r 12 , (83b) 



the spectrum [Eq. (fSTj) ] as well as the intensity correla- 
tion [Eq. ([55])] associated to Eq. fl70j! [and Eq. fl?!}] are 
indistinguishable from those shown in Fig. 5. In spite of 
these similarities, both cases lead to very different photon 
counting processes. 

In Fig. 6, the photon counting statistics is character- 
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ized through the stationary Mandel factor Q st [see Ap- 
pendix C] as a function of the system-laser detuning, i.e., 
S = lol—oja [Eq. (|26|) ]. While the number of emitted pho- 
tons (the stationary intensity) is inversely proportional 
to 6, in the limit of large detuning the Mandel factor 
may take arbitrary values. In fact, this object [Eq. ([68]) ] 
measures the (normalized) fluctuations in the number of 
counted photons around its mean value. For the light as- 
sisted dynamic [Eq. (f76|) ] (solid line), the photon statistic 
becomes super-Poissonian, i.e., Q s t > 0, with the asymp- 
totic value 



General environment fluctuations 



Um q t = 2 7i27 2 i [(7i + 72i) ~ (7 2 + 7i 2 )] 2 
s ->°° St (7i2 +72i) 2 (7i7i2 +7 2 72i + 27i 2 72i) 



> 0. 

(84) 

Notice that this expression does not depend on the ex- 
ternal laser excitation. This result is expected because 
the condition ri <gc <5 is satisfied. Instead, the photon 
counting statistic for the self fluctuating environment 
[Eq. ([70]) ] (dotted line) becomes Poissonian for larger de- 
tuning, lima^oo Q s t — 0. On the other hand, both dy- 
namics are characterized by the same Mandel factor when 
the system-laser detuning is small. The inset shows the 
line shape Eq. ([57]) associated to the light assisted pro- 
cess. 

The super-Poissonian character [Q s t > 0] of the pho- 
ton counting process in the limit of increasing system- 
laser detuning has a clear origin. For the light assisted 
process, all transition rates Tr>r [Eq. l)82p] vanishes in 
the limit S — > oo. Instead, for the self- fluctuating envi- 
ronment the transition rates does not vanish in the same 
limit. In fact, they are independent of the properties of 
the external laser excitation, Eq. (fTTj) . The vanishing of 
lima^oo Qst can be recover when the asymptotic value of 
one of the transition rates [Eq. ([B"2")) ] does not depend on 
the system-laser detuning. This property can be achieved 
by assuming the following scaling (dashed line in Fig. 6) 



712 -> 7i2 + 7i 2 (l<5IA>), n -> si+n(\s\/8 ) 1/2 , 

(85) 

where So, f2, and 7 12 are arbitrary constants. Both 
7 12 and f2 increase as S increases. Consequently, the 
transition rate Ti2 does not vanish by increasing <5. In 
fact, lim5_ +oo ri2 ~ 7i2^ 2 /(7i2 + 4<5q). Furthermore, as 
72 ^ 7i2i "given" that the bath is in the configurational 
state 1 2) , most of the photon emissions are attempted by 
the transition |2) — > |1) , implying 7 2 — 0. Therefore, with 
the scaling Eq. (|85p . the dynamic can be mapped with a 
triplet blinking modelling [J |46|, [52| where the configura- 
tional state 1 2) defines a dark state incoherently coupled 
to the system. On the other hand, by assuming valid the 
condition 7 21 <C 7i, from Eqs. ([75)) and (|52")l it follows 



h ~ 7l ft 2 /(4<5o|<5|), and T 21 ~ j 21 tt 2 /{45 \5\), which in 
turn implies that the bath state |1) can be associated to 
the bright intensity states. Finally, from Eq. l)84p. the 
scaling Eq. ([83)1 implies lim^oo Q st ~ 2721/71 ^ 1- 



Having understood the derivation of Eq. ([19)) and their 
associated physical processes, one can quickly write down 
master equations that introduce more general environ- 
ment fluctuations. The evolution of the auxiliary states 
is written as 



dp R {t) 
dt 



&U4?)[Pfl(*)]- 



(86) 



The superoperators effi and £^ define the unitary 
and irreversible evolution for each configurational bath 
state respectively. Both kind of contributions have been 
analyzed in the previous examples. Nevertheless, notice 
that extra pure dispersive contributions (phase destroy- 
ing process) may also be considered. The non diagonal 
superoperator C jjr introduces the more general environ- 
ment fluctuations. It can be written as 



J cf [Pr 



] = -E ^{^w+- 



Vrr' Ap R ,A\ 



R' 
R'^R 



E 

R' 
R'^R 

(87) 

defines the characteristic non- 



A is an arbitrary system 



where the matrix rj R i R 
diagonal rates of the problem, 
operator. 

Only when A = I, where the I is the identity oper- 
ator, the environment fluctuations, defined by the evo- 
lution of Pr{€), do not depend on the state of the sys- 
tem, recovering Eq. ([70)) . When A = a = \a) {b\ , the 
evolution Eq. ([75)) is recovered, which is associated to 
bath fluctuations that are attempted by the system tran- 
sition \b) — > I a) , i.e., by a photon emission. When 
A = er T = \b) {a\ , the configurational transitions are at- 
tempted by a system transition from the lower to the 
upper state, \a) — ► \b) , i.e., the system absorbs a photon 
from the background electromagnetic field. This kind of 
contributions, which also arises from the microscopic in- 
teraction Eq. were discarded because at room tem- 
peratures the optical thermal excitations are almost null. 

For A = \b) (b\ , Eq. ([87)) describes configurational tran- 
sitions that can only happen when the fluorophore is in 
the upper state. These processes can be microscopically 
described by replacing the interaction Hamiltonian Huw, 
Eq. (fT5)) . by \b) (b\ <g> H uw . They induce light assisted 
processes similar to those described by Eq. ([76")) . In par- 
ticular, when 6 — * 00, the stationary photon counting 
statistics is also super-Poissonian. Similarly, A = \a) (a\ , 
describes configurational transitions that only happen 
when the system is in the ground state. In both cases, 
the dynamic induced by Eq. (J87j introduces a dephasing 
mechanism that affects the fluorophore yield. 

The formalism does not forbid the situation in which 
Vri-r = 0. This case was (partially) addressed in Refs. 
[5 lL |52| . where the underlying dynamic was defined in 
terms of a structured reservoir whose influence can be ap- 
proximated by a direct sum of Markovian sub-reservoirs 
(generalized Born-Markov approximation) [4^, |4^| . Most 
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of the expressions for the scattered field observables ob- 
tained previously remain valid. Nevertheless, some prop- 
erties of the radiation pattern may depend on which way 
the system-reservoir fluctuations are measured (59[. 

VI. SUMMARY AND CONCLUSIONS 

The main theoretical problem of SMS is to relate the 
radiation of a single fluorophore system with the under- 
lying dynamical fluctuations of its local environment. In 
this paper, we tackled this problem from a quantum- 
electrodynamic treatment based on an open quantum 
system approach. 

The central ingredient of our formalism is the de- 
scription of the nano-environment. Instead of consid- 
ering the microscopic dynamics associated to each spe- 
cific situation, after noting that the Hilbcrt space struc- 
ture of the reservoir can not be resolved beyond the 
experimental resolution, it was approximated by a set 
of configurational macrostates, each one taking in ac- 
count all microscopic bath configurations that lead to the 
same system dynamic. This coarse grained representa- 
tion allows to define an effective microscopic description, 
where the fluorophore and background electromagnetic 
field Hamiltonians, as well as their mutual interaction, 
are parametrized by the configurational states. From 
the effective microscopic Hamiltonian, it was possible to 
achieve the central goals of this paper. 

By tracing out the electromagnetic field and the con- 
figurational degrees of freedom, we derived the Lindblad 
rate evolution Eq. (|19p . which encode the statistical be- 
havior of the fluorophore and the reservoir fluctuations. 

In contrast with previous approaches, the scattered 
electromagnetic field was characterized from a full 
quantum-electrodynamic description. The electric field 
operator was written in terms of conditional system op- 
erators, Eq. (f3l>)) . After appealing to a quantum regres- 
sion theorem, they allow to expressing the field correla- 
tions in terms of the system density matrix propagator, 
Eqs. P5")l and (|4"6l) . High order correlations can also be 
expressed in a similar way. These results allow to get 
simple manageable expressions for observables like the 
spectrum, Eq. (|48|) . and the intensity- intensity correla- 
tion, Eq. ([55]). 

An extended Mandel formula, Eq. (fSl)|) , which includes 
an additional sum over all configurational transitions, de- 
scribe the photon counting statistics. A simpler charac- 
terization based on the generating operator Eq. ([50)) was 
presented. The line shape and stationary Mandel factor 
can be obtained straightforwardly after solving its evolu- 
tion, Eq. (|6Tj) , in the Laplace domain. 



The examples worked out in the manuscript were clas- 
sified according to the properties of the configurational 
fluctuations. Processes like spectral diffusion and lifetime 
fluctuations correspond to bath fluctuations that do not 
depend on the system state while light assisted processes 
are recovered in the opposite case, i.e., when the environ- 
mental fluctuations become entangled with the system 
dynamic. From a microscopic point of view, the former 
case arises in local environments characterized by dense 
manifold of states while the last one arises from reservoirs 
defined by few degrees of freedom. A recipe for studying 
arbitrary bath fluctuations is given by Eq. (|57|) . 

The approach allowed to study observables and phe- 
nomena that cannot be easily obtained from formalisms 
that do not take explicitly into account the quantum na- 
ture of the scattered electromagnetic field. We showed 
that the incoherent optical spectra may develops nar- 
row peaks whose width allow to read the value of some 
of the characteristic rates that govern the reservoir fluc- 
tuations. A convergence to a super-Poissonian photon 
counting statistic (in the limit of large system-laser de- 
tuning) in light assisted processes was analyzed in detail. 

The open quantum system approach provides a solid 
basis for modeling a broad class of SMS experiments. 
The effective microscopic dynamic and the quantum- 
electrodynamic treatment give the possibility of explor- 
ing many open problems. Since a density matrix de- 
scription is available, a general formulation of a quantum 
jump approach should provide insight into the properties 
of the photon-to-photon emission process. The incidence 
of non-stationary phenomena and quantum effects in the 
configurational space are also interesting situations that 
can be dealt with the ideas introduced in this paper. 
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APPENDIX A: GENERALIZED OPTICAL 
BLOCH EQUATION 

The notation of the matrix elements Eq. (f52|) can be ex- 
tended trivially to each contribution Gn(t 7 s) in Eq. (JSpJ- 
The evolution of their matrix elements, i.e., the general- 
ized optical Bloch equation, from Eq. (f6Tj) . read 
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U R (t,s) = S R V R (t,s)- -l R + ^(t> R - R U R {t,s) + ^(j) RR ,U R '{t,s), 

w R (t,s) = n R v R (t,s)-^ R + s 7R )[w R (t,s) + y R {t,s)] 

- 1 E ( f > s ) + (*> s )i - E ^H'fl w fl(*: s ) + E ^rr' w «' (*> s ) 



s) - n^w^t, s) - [-7* + E ^K' J V ^*' s ) + E ^RR> V R' (*> s 



(Ala) 
(Alb) 
(Ale) 



3>k(M) 



5 E ^ ( f > s ) + y R- (*> s )i - E ^H'fl^(*» s ) + E ^yfl' (*, «), 



(Aid) 



r 



where 5 R = (ljl — u>a)—8u>a Pi- 0], and 7 fl being The contributions (d/ds)G R (t, s)| s=1 , can also be ob- 
defined by Eqs. iJUJ and respectively. 



APPENDIX B: INTENSITY-INTENSITY 
CORRELATION IN THE GENERATING 
FUNCTION APPROACH 

In this Appendix, the intensity-intensity correlation 
Eq. (|53j) is obtained from the generating function ap- 
proach Eq. (jnj). These calculations shows the consis- 
tency between the generating function approach, the 
Mandel formula Eq. (|56[) . and the results of Sec. III. 

The function y(t,s), Eq. (|62d[) . from the definition 
Eq. ([6"0|) , can be written as 

y(t,s) = hrr s [G(t,a)] = \^2 R Tr s [G a (t,a)]. (Bl) 

Then, the evolution Eq. (fuTj) allow us to calculate 
the time derivative of the factorial moments N^ n \i), 
Eq. ([64)) . For n = 1, we get 



-iV W(t) = £ 7«Tr 5 [a Wi)] = ]T 7* (&| 1 6) 



R 



(B2) 



Since iVW(i) = N(t) [Eq. (|66a|) ]. the line shape Eq. J67] 
can alternatively be written as 



/K) = E7*<&Ip£ |&) 



(B3) 



where = hmt—Kx, p R (t). This expression recovers 
Eq. (|55|) . showing that the line shape is proportional to 
the initial value of the first order correlation, Eq. ([37)) . 
The same expression follows straightforwardly from the 
Mandel formula, Eq. ([56")) . By using the same procedure, 
the evolution of the second factorial moment reads 



|A( 2 > W =2]T 7fl Tr 5 







cr T er —G R (t, s) 



(B4) 



tained from Eq. (JHI 
d d 



dt ds 



G R (t,s) 



= ^2£rr' -~-G R >(t,s) +a R (t), 
r> ° s »=i 

(B5) 

where is defined by Eq. (|4"4")) and the inhomogeneous 
term a^;(t) reads 



After a formal integration of Eq. (|B5|) . from Eq. (|B4 
N^ 2 '(t) can be written as 



N<- 2 \t) 



dt' 



dt"F{t' -t",t"), 



(B7) 



where 



F(r, = E 7fl^s{ffV(e T ^)jiB, [a* (<)]}. (B8) 



Notice that this expression, with the exception of the 
spacial and angular dependences, recovers Eq. (|46p . Now, 
from the Mandel formula Eq. ([56]) . after a standard set 
of calculations steps [f| B H9] > it follows the relation 



N^(t) 



dt' 



dt" :I(t'- t")I(t") 



(B9) 



Therefore, the intensity-intensity correlation reads 



9 {2) {r) 



lim t 



r(r,t) 



(BIO) 



which from Eq. (|B8[) recovers Eq. ([53")) . These results 
show the internal consistency of the developed approach. 



17 



APPENDIX C: STATIONARY MANDEL FACTOR 

In this Appendix, the stationary Mandel factor 
Eq. (55]) is obtained from the function y(t, s) [Eq. ([52])] 
in the Laplace domain. From Eq. ([55]) , the Mandel factor 
can be rewritten as 



Q(t) = 



y"(t, l) - 2[y f (t, i)} 2 



(Cl) 



where the accent denotes derivation with respect to s. In 
the long time limit, the derivatives of y(t, s) behave as 



y'(t,l) w a + bt, 
y"(t,l) w C + At + Bt 2 . 



(C2a) 
(C2b) 



The constant a is zero only when the system begins in its 
stationary state. From Eqs. (|C1|) and (|C2[) the stationary 
Mandel factor, under the condition B = 2b 2 , reads 



Qst = -r ~ 4a. 
o 



(C3) 



The condition P = 2b 2 guarantees the cancellation of 
terms cuadratic in time. On the other hand, from the 



definition Eq. ([67]) . the line shape reads 
J(wi) = 2b. 



(C4) 



The coefficients associated to the asymptotic time behav- 
iors, Eq. (|C2[) . can be obtained in the Laplace domain. 
In the limit of u — > 0, it is always possible to write 



y'(u,i) » 

implying the relations 
6 = P' 

25 = £ 

P 



p + qit + 0[m 2 ] 
Pu 2 + Qu 3 + 0[u 4 } ' 

p + qu + 0[u 2 ] 
Pit 3 +Qu 4 + 0[u 5 }' 



q pQ 
p P 2 ' 



(C5) 
(C6) 

(C7) 



A=l-B. (C8) 
p pi y ' 



Then, the line shape and stationary Mandel factor can 
be obtained in an exact way after solving the generalized 
optical Bloch equation [Eq. (|A1[) ] in the Laplace domain. 
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